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We study different solid phases of 4 He, by means of Path Integral Monte Carlo simulations based 
on a recently developed worm algorithm. Our study includes simulations that start off from a high-T 
gas phase, which is then "quenched" down to T=0.2 K. The low-T 1 properties of the system crucially 
depend on the initial state. While an ideal hep crystal is a clear-cut insulator, the disordered system 
freezes into a superglass, i.e., a metastable amorphous solid featuring off-diagonal long-range order 
and superfluidity. 
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The remarkable observation by Kim and Chan of a 
non-classical moment of inertia in solid 4 He Q has gen- 
erated a new wave of interest in the possible superfiuid 
phase of a solid. Supersolidity of 4 He is still controver- 
sial, both at the experimental and theoretical levels. Two 
of us have recently proven that, irrespective of its mi- 
croscopic structure, any supersolid crystal should con- 
tain gapless vacancies and/or interstitials In other 
words, any continuous-space supersolid is generically in- 
commensurate (i.e., the number of atoms per unit cell is 
not an integer) and squeezeable, i.e., by applying pressure 
it should be possible to squeeze matter from a container 
with supersolid into a buffer volume containing the same 
supersolid. However, this very experiment has yielded a 
negative result for solid 4 He y|- 

A wealth of numerical studies clearly indicate that 4 He 
is a commensurate (thus insulating) crystal. The finite 
activation energy of a vacancy computed numerically is 
large, ~ 15 K, and claimed consistent with the experi- 
mental observations Q. The activation energy of an in- 
terstitial, ~ 50 K H, is significantly larger than that of a 
vacancy. A simulation study of exchanges in an ideal hep 
crystal []|, yielded indirect evidence that the system is 
not superfiuid. In sharp contrast, the variational (T=0) 
calculation of Ref. claims a finite condensate fraction 
in the commensurate 4 He crystal. Thus, additional in- 
vestigation is warranted. 

The experiment of Kim and Chan itself has revealed 
a number of facts pointing to a strongly inhomogeneous 
scenario of superfluidity, chiefly the contaminating effect 
of a small concentration of 3 He, and non-XY behavior of 
the superfiuid density at the critical temperature. The 
need of exploring inhomogeneous (metastable) scenar- 
ios of supersolidity, dictated both by theory and experi- 
ments, has already resulted in some relevant theoretical 
developments, revealing superfiuid interfaces in a lattice 
solid |jj and a superfiuid layer at the boundary between 
the 4 He crystal and a disordered substrate jjj]. 

The numerical observation of a metastable disordered 



supersolid, (a superglass phase of 4 He) is reported in the 
present Letter. To be specific in the definition, by glass 
we mean a spatially disordered (metastable) phase, in- 
distinguishable from a solid on a time scale much 
shorter than the typical relaxation time, i re i, which in 
turn should be dramatically longer than the inverse De- 
bye frequency, w^ 1 ■ Superglass is the term that we use for 
such a phase, if it also displays superfluidity. Note that 
our definition of glass does not address the behavior of the 
system at time scales t > ireh whereupon it may undergo 
structural relaxation into the polycrystalline sample, or 
simply behave as a very viscous liquid. 

In Ref. 0| , the idea of glassification of overpressur- 
ized liquid 4 He was put forward, in order to explain a 
striking experimental outcome, i.e., the absence of bulk 
solid nucleation under fast (about 1 /is) acoustic wave 
compression pulses, up to pressures as high as ~ 160 bar. 
The authors conjectured that the glassy phase is normal 
(though the experiment was done at T = 0.05 K and the 
adiabatic heating was estimated to be below 0.1 K); the 
absence of superflow towards the nucleation center would 
dramatically suppress the rate of growth of the crystal. 
Conceptually, the finding of the present Letter is differ- 
ent, but we believe relevant to the interpretation of the 
experiment of Ref. [1(1 ] . Jamming of structural relaxation 
does not per se exclude superfluidity. Crystallization is 
suppressed by the mere fact that the normal component 
forms a glassy solid, implying that further evolution to- 
wards a lower-energy polycrystal structure necessarily in- 
volves a chain of exponentially rare quantum-tunnelling 
or thermoactivation events, rather than a rapid growth 
of the supercritical nucleus. Indeed, the boundary be- 
tween the perfect-crystal nucleus and the superglass, is 
a solid-solid interface which realizes a pronounced local 
energy minimum. Its evolution should therefore imply 
either quantum tunnelling (in the T — > limit), or ther- 
moactivation. 

Our study is based on accurate Path Integral Monte 
Carlo (PIMC) simulations of condensed 4 He, making 
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FIG. 1: (Color online). Pair correlation function of the ideal 
4 He hep solid and superglass at the near-melting (n=0.0292 
A -3 , lower panel) and higher (n=0.0359 A -3 , upper panel) 
densities. 



use of a recently developed worm algorithm 11]. This 
method allows for efficient sampling and accurate deter- 
mination of the single-particle Green function and su- 
pcrfluid density, for systems comprising a relatively large 
number TV of particles (of the order of several thousand). 
Specifically, we address the following two issues: (i) Is 
it possible to obtain definitive first-principle theoretical 
evidence that an ideal hep 4 He crystal is an insulator ? 
(ii) What happens to a sample of liquid 4 He quenched 
through the first-order liquid-solid phase transition? 

We consider a system of TV 4 He atoms (7V=216 and 
800), at a temperature 0.2 K < T < 1 K, and at the 
two densities n=0.0292 (0.0359) A -3 , corresponding to 
an ideal hep 4 He crystal at a pressure of approximately 
32 (155) bars The sample cell geometry with peri- 
odic boundary conditions is designed to fit an ideal hep 
crystal. We use the standard microscopic model of 4 Hc, 
based on the Aziz pair potential [l3| . 

In Figs. Qand El we show data for the pair correlation 
function g(r) and the single-particle density matrix n(r). 
For both the near-melting density of n=0.0292 A~ 3 and 
the higher density of n=0.0359 A -3 we study two sam- 
ples, differing in one respect only, namely their initial 
configurations before equilibration. 

The single-particle density matrix is defined as 
n(r, r') = (V>'(r) tp(r')) where ^(r) is the particle an- 
nihilation operator, p(r) = ^(r)-0(r) is the local 4 He 
density operator, and (...) stands for thermal average. It 
is customary to display the spherically averaged function 



n(r) 



— [ dQ f d 3 r' n(r',r' + r) 
AttV J J 



(1) 



where V is the volume of the system. This is the quantity 
shown in Fig. [21 

When the simulation is started from an initial config- 
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FIG. 2: (Color online). Single-particle density matrix n(r) 
for ideal 4 He hep crystal (filled symbols) and superglass (open 
symbols) at the near-melting density n,=0.0292 A~ 3 (squares) 
and high-density n=0.0359 A -3 (circles). Solid lines through 
filled symbols represent exponential decay. 



uration corresponding to an ideal hep crystal, we con- 
sistently find an exponential decay of n(r) at large dis- 
tances, with short-range oscillations due to coordination- 
sphere effects. We observe no change in the results be- 
tween the temperatures of 0.2 and 1 K, to indicate that 
those shown in Fig. are essentially ground state esti- 
mates. This result provides a robust confirmation that 
an ideal hep crystal is not a Bose condensate (superfluid) . 
This conclusion is consistent with the theoretical expec- 
tation that a crystal with finite activation energies for 
vacancies and intcrstitials will not display superfluidity 
0, and is in agreement with arguments based on the 
statistics of exchange cycles observed in the same system 
0| . The results and conclusions of Ref. |(| appear there- 
fore to be erroneous, possibly artifacts of the variational 
approach. 

At present, there is no clear understanding of what 
crystalline defects dominate in the experimental samples 
of Ref.[lJ. It is not known whether individual dislocations, 
dislocation sheets and networks, or grain boundaries in 
bulk 4 He may underlie the experimentally observed su- 
perfluid response, though model simulations of domain 
walls in quantum solids hint at such possibilities 7]. But 
regardless of their nature, in the absence of crystalline 
defects no theoretical interpretation seems viable of the 
experiments reported in Ref. [l[ in terms of superfluid 
response. 

In order to investigate scenarios of broken transla- 
tion invariance not involving a perfect crystal (though 
perhaps not directly related to the experiment of Ref. 
0J), we designed a simulation protocol aimed at mimick- 
ing a "quenching" experiment (namely one in which liq- 
uid Helium is suddenly, rapidly cooled) obviously mak- 
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ing allowance for the important differences between the 
imaginary-time PIMC dynamics and the real-time dy- 
namics of actual physical systems (see discussion below) . 
Starting from an initial configuration characteristic of a 
high-T gas phase, we "quench" the system down to the 
temperature T = 0.2 K by > 10 4 PIMC sweeps. One 
"sweep" is defined as the number of accepted updates 
sufficient to sample the entire Path Integral configura- 
tion. We then run the simulation long enough to achieve 
stability of statistical averages for structural properties 
and for the single-particle density matrix. 

The phase that emerges from disorder resembles the 
hep crystal only at short interatomic distances (of about 
two-three coordination spheres), with no diagonal long- 
range order (see Fig. Moreover, this phase has a 
well-developed off-diagonal long-range order, with a con- 
densate fraction n D = 0.5 % (see Fig. 01, and a surpris- 
ingly large superfluid fraction p s = 0.6(1) at n=0.0292 
A" 3 and p s = 0.07(2) at n = 0.0359 A~ 3 . Though the 
superfluid fraction is strongly suppressed with pressure 
the condensate fraction is reduced by merely 50%. 

The important observation is that 4 He can remain in 
the metastable superfluid state at solid-state densities, 
even at fairly high pressure. This observation is consis- 
tent with a previous study |l4j . predicting overpressur- 
ized liquid He to remain superfluid (at T=0) to arbi- 
trarily high density. The nature of the superfluid phase 
will generally depend on pressure, temperature, and the 
experimental time scale. For example, one may expect 
that the lower density finite-temperature phase should 
be just a superfluid, but with a rather viscous normal 
component. On the other hand, with increasing density 
such a normal component may evolve into a glass, with 
a diverging (i.e., unobservably large) viscosity. 

In order to study whether and how the system breaks 
translation symmetry, we calculate the condensate wave 
function 4>(r). The worm algorithm offers direct access to 
the one-body density matrix which in the presence of off- 
diagonal long-range order factorizes at large separation 
|r-r'| 

(^(r)^(r')) -> 0(r)0(r') . 

In Fig. we show two-dimensional xy-maps of the 
condensate wave function cj){x,y,z) at n = 0.0359 A~ 3 , 
for ten (equally spaced) slices along the z-axis. The data 
shown in Fig. [3] represent long simulation-time averages, 
not instantaneous snapshots. Aside from the obvious ob- 
servation that the system has manifestly broken trans- 
lational invariance, the results suggest no obvious inter- 
pretation of the disordered pattern for <f>(r) in terms of 
dislocations or grain boundaries . Correlations barely 
extend over two slices in Fig. [21 The conclusion that we 
draw from Figs. 11131 is that 4 He forms a superglass. 

Naturally, the results shown in Fig. [3 are influenced 
by a particular gas-like initial condition; another ini- 
tial condition would produce a different result for 4>(r). 



Nonetheless, the fact that superfluidity and off-diagonal 
long-range order appear for just one such random ini- 
tial condition, strongly suggests that these are genuine 
physical properties of the metastable phase. 

Despite the above-mentioned fact that the MC dy- 
namics is quite different from the real-time dynamics of 
helium (in simulations heat dissipates locally) it is pos- 
sible to make semiquantitative arguments with regards 
to the stability of the superglass phase. Its stability 
on timescales several orders of magnitude longer than 
tap 1 ~ 3~ 13 s, is guaranteed by observing no changes 
in the superglass properties over 10 4 MC sweeps. [The 
most appropriate physical interpretation of one sweep for 
the conventional PIMC scheme, is the time scale cor- 
responding to the zero-point motion of atoms, whereby 
all particles have a chance to sample their optimal po- 
sitions locally.] The extra advantage of using the worm 
algorithm is that pair-wise exchange, or tunneling of two 
particles, is sampled at the same rate as zero-point vibra- 
tions, while the rate J of exchange processes in the solid 
state of 4 He is about five orders of magnitude slower (e.g., 
measured values of tunneling for 3 He in 4 He are of order 
J ~ 10 ps -1 0|). It seems then plausible to assume 
that the metastability of the superglass phase extends 
up to 10 4 J -1 ~ 1 ms. If multiparticle tunneling (at low- 
temperature) events are required to reach the genuine 
equilibrium, then the actual degree of metastability has 
no obvious upper limit, and can easily exceed the exper- 
imental timescale. 

Summarizing, we have provided theoretical evidence 
that 4 He features a new metastable phase, a superfluid 
glass. This observation naturally suggests that other, 
more "regular" types of solid disorder, such as grain 
boundaries and dislocations, may also possess superfluid 
properties. We foresee further theoretical studies in the 
following directions: 

1. Determination of the "phase diagram" of the glassy 
phase, i.e., of its domain of metastability (e.g., in 
the (n, T) plane) and the line separating superglass 
from normal glass. The other important issue is to 
quantify the crossover line separating this novel, 
superglass phase from the superfluid: the super- 
glass phase is characterized by the low-temperature 
plateau p s {T) — > p < 1, as in a "dirty" conven- 
tional superfluid. 

2. Explore alternative possibilities of metastable su- 
persolids: (i) A regular crystal doped with vacan- 
cies (or even interstitials) . (ii) Superfluid grain 
boundaries and/or dislocations, (iii) The yet elu- 
sive superfluid phase of condensed para-hydrogen. 

The superglass phase is not directly relevant to the in- 
terpretation of Kim and Chan's experiment, since the 
MC temperature quench is much more rapid then in 
the experiment, leading to more disordered samples and 
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FIG. 3: Condensate wavefunction at n=0.0359 A -3 (repre- 
sented by the density of points) and T=0.2 K, obtained by 
making ten slices of the system along the ^-direction and pro- 
jecting them on the a;i/-plane (slices are ordered from left bot- 
tom to left top and then from right bottom to right top). 



much larger superfluid fraction. However, a related sce- 
nario may be appropriate, namely that of a generalized 
superfluid-grain-boundary [tJ, which may may include 
a foam-shaped superglass network interpenetrating the 
polycrystal. 
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Note added. In an independent study [T^ |. Clark and 
Ceperley calculated the single-particle density matrix of 
the ideal hep 4 He crystal at the melting density n — 
0.0287A" 3 . The data of Ref. [13 are consistent with 
ours for n = 0.0292A" 3 . 
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